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Abstract 

We describe the origins of recurrence relations between field theory amplitudes in terms of the 
construction of Feynman diagrams. In application we derive recurrence relations for the ampli¬ 
tudes of QED which hold to all loop orders and for all combinations of external particles. These 
results may also be derived from the Schwinger-Dyson equations. 


1 Introduction 

There is currently a wealth of interest in the perturbation series of Yang-Mills theory expressed 
in terms of maximal helicity violating, or MHV, amplitudes. These amplitudes may be used as a 
set of vertices with which to generate an alternative perturbation expansion of the theory. 

However, until very recently there was a tree level proof that these results hold generally 
HI, with numerous examples at low loop order g - J. It has now been shown how to write 
the Yang-Mills action in terms of variables which makes the MHV structure manifest The 
quadratic part of the action is the inverse of a scalar propagator which connects positive and 
negative helicity fields while the infinite number of tree level MHV amplitudes appear as an 
infinite number of vertices in the action. This is an exciting result since it allows perturbative 
Yang-Mills to be re-developed in terms of Feynman diagrams constructed from MHV amplitudes. 
This approach describes amplitudes through recurrence relations, a computationally more efficient 
approach than the usual Feynman diagram expansion. Examples of similar recursion results have 
been found to hold in QED and gravity El ED, though only at low orders. 

There is, then, considerable interest in alternative perturbative structures in quantum Held 
theory. In this paper we will show that such recursion relations may be derived from properties 
of correlation functions. We will construct recursion relations for QED amplitudes which hold 
for all external particles and to all orders in the loop expansion. 

This paper is organised as follows. We begin with an example, using scalar fields, which 
illustrates our method. The example is based loosely on arguments from the sewing of string 
world sheets and suggests that the object to be calculated is the expectation value of the free 
action. In section 0 we evaluate this in terms of functional operators and use power counting 
arguments to generate recurrence relations between correlation functions to all loop orders. We 
explain how this leads directly to recurrence relations between scattering amplitudes. In section 
Q] we extend our arguments to QED, deriving recurrence relations for the scattering of photons 
and fermions. Finally in section |S] we give our conclusions, summarise an alternative derivation 
from the Schwinger-Dyson equations and discuss the extension to other theories. 
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2 An example 

Consider the four particle scattering amplitude in rj> 3 theory, given by the LSZ formula 

(Pi,P 2 |P3,P4> = f dxi...d® 4 e lpi X1 e W2 ' X2 e ~ %P3 ' X3 e ~ lpiXi '^\i{d 2 x . +m 2 )( 1...4) (1) 

J j=i 

where the pi are on shell. Our metric is mostly minus. We denote an integration over all D + 1 
dimensional spacetime by dx and we have abbreviated the connected correlation functions to 

( 1 . . . n ) = { 0 | T[4>(x l) . . . 4>{x„)\ | 0 ) connected ■ 

In scalar </> 3 theory the tree level contribution to the correlation function above is 

(1... 4 }o = —A 2 fdzdz G xlz iG X2Z iG zz iG ZX3 G ZXi (2) 

£7 

where the subscript zero indicates tree level, G xy is the propagator between points x p and y p 
and the sum over a is over inequivalent permutations of the external points, in this case giving 
three possible diagrams each with a factor of —A 2 , below. 





The key observation is that these diagrams may be split into three field diagrams be severing the 
internal line. This can be performed explicitly using the particle equivalent of Carlip’s sewing 
method for string worldsheets 0. In the string case two worldsheets are sewn together at a 
boundary by integrating over shared (co-ordinate and ghost) boundary data but with an insertion 
of the string Hamiltonian acting on one of the boundaries. This is necessary to correctly cover 
the moduli space of the worldsheet being formed. 

The same consideration must be applied to particle theories, albeit at a basic level. We 
will be interested in cutting and joining propagators. A careful derivation of the propagator 
which maintains reparametrisation invariance, see 0, reveals a modular parameter which is the 
intrinsic length of the worldline. This is usually not manifest since the propagator is commonly 
represented by its Fourier transform. The propagator here is G xy , the inverse of i5(x — y)(d x + m 2 ) 
(the Gaussian operator in the partition function with an appropriate ie prescription). For this 
object Carlip’s sewing becomes the almost trivial 

G xy = J dz G xz i(d 2 + m 2 )G zy . (3) 

Returning to our example we apply 0 twice to the internal line, separating the four held diagrams 
into two three held diagrams. The correlation function can then be written 


<l-..4>o = £ J 


dadfe 


i(dl + m 2 )( 1, 2, a)o 


Gab 


i(d 2 +m 2 )(3,4, b) 0 


(4) 


where the sum over a takes in the symmetry of G a b and therefore the equivalence of the two 
terms in large brackets. This may be checked using the explicit form of the amplitude in ©. 
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The insertions of the Hamiltonian are of course the operators appearing in the LSZ reduction 
formula. If we insert our result into 0 and write G a b as its Fourier transform, 


Gab — ^ 


/ 


d D+1 k e ~ ik G~b) 


(27r) D+1 k 2 — \ 


which provides the necessary exponential factors, then the four particle scattering amplitude 
becomes 


r dfc i 

<Pl,P2|P3,P4>0 = X, J ( 27r )D+l (Pl.P2| fc )0^J—^(fc|P3,P4) 0 (5) 

where the sum is now over inequivalent permutations of the external momenta (including k). 
Since each S matrix element is proportional to a momentum conserving delta function we may 
perform the integral over leaving 


A(j)l,P2 |P3,P4)o = A{pi,p 2 | k ) c 


k 2 — m 2 


A{k |p 3 ,P4)o 


( 6 ) 


where k conserves momentum, k = pi + P 2 = P 3 + P 4 - We have expressed the four particle 
scattering amplitude in terms of three particle scattering amplitudes connected by a propagator 
(a structure familiar from, but not the same as, the MHV rules). In the next section we will 
derive a general form of this result which holds to all loop orders. 


3 Scalar fields 


Before discussing QED we will give the general result of applying Carlip’s method to a scalar Held 
theory with interaction \(j> N /N\ We would like to show that all connected n-point correlation 
functions are related to n + 2-point correlation functions via the use of Let the partition 
function for the theory be Z[J] = where W is the generator of connected correlation 

functions. We are interested in proving something like 


W[J] ~ Jdxdy (dl + m 2 )^^G xy (d 2 y + m 2 ) ^^ 


h 


<5W[J],„ 2 2 JW[J] 


SJ(x) 


{d x + m ) 


SJ(x) 


(7) 


where the two additional fields on the right hand side become derivatives with respect to the 
source J. The second line in the above follows from an integration by parts and points us at the 
object we will investigate, namely the expectation value of the free action So in the presence of 
a source, 

k So \~ ih jj\ z [ J ] = J T>< ^ I So ^ exp + \ J ( 8 ) 

Introducing a parameter £ we can write 0 as 




/ D0 exp f 

C=o J V 


«(i + 0 

h 


So[0] + l jS N [<t>] + 


(9) 


The free scalar field action is 

5o[0] = i J dx d4>.d(j> — m 2 4> 

so the left hand side of 0 is 



i{dl + m 2 ) 


8 2 W 

5J(x) 2 


+ 


5W 

8J(x) 


i(d 2 x + m 2 ) 


5W 
8J(x)' 


( 10 ) 
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To evaluate the right hand side of ©, observe that W is the sum of connected diagrams con¬ 
structed from V vertices, n external lines and I internal lines, with l loops. We expand W[J\ as 

1 _°°. C n 71 

W[J] m J / dxi . . . dx n —J{x l) . . . J(x n )( 1 . . . tl). (11) 

n=0 ** n ’ 

Each of the correlation functions has a loop expansion in powers of h, 


OO 

{l...n) = J2 hl ( 1 ---n)i (12) 

1=0 

with l = 0 being the tree level pieces etc. The only change to the generating functional W which 
results from the inclusion of the £ term in the exponential is that each propagator is accompanied 
by a factor (1 + £) _1 . Therefore each Feynman diagram picks up a factor of (1 + C) - " -7 - Using 
the Euler formula and the loop counting formula, 

n + 2I = NV, I-V = l- 1, (13) 


we deduce that each set of connected Feynman diagrams which make the correlation functions in 
the expansion of W now carries a factor 

(1 + £) -n-J l n ’d where I(n,l) \= (n + N(l — 1))/(N — 2). (14) 


Returning to the derivative with respect to ( brings down this new sum of connected diagrams 
from the exponential. Setting £ to zero leaves W with — n — I(n, l) multiplying the n-field 1-loop 
correlation function, all multiplied by the usual partition function. Equating both sides of we 
find, 


^ Jdx i(dl + m 2 ) 


SJ(x) 2 5J(x) 


i(d% + m 2 ) 


8J(x) 


°°^ r jn °°^ 

J2 / dxi... dx n —J(x i)... J(x n ) h\n + I(n, l))( 1 ...n)i. 

-k,— n J ;—n 


(15) 


Equating the coefficients of J n h l 1 gives 

(n + I(n, l))( 1 .. .n); = i Jdx i(dl + m 2 )( 1... n,x,x)i 


+ 2 S y m 


nip 

EE/ 

P =0 l'=0 J 


A /, , i{dl + TO 2 ) , , 

dx(l...p,x)v— - —(p + 1 ...n,x)i_v 

p\(n — p)\ 


(16) 


where the symmetrisation is over all possible distributions of x\ to x n in the pair of correlation 
functions. This is our result. It describes recurrence relations, in the number of particles and the 
loop order, between the correlation functions of scalar field theory which hold to all loop orders. 


To be thorough, it can be seen that when the left hand side is (1 . . ,n)i the right hand side 
contains terms multiplying correlation functions of n fields and n + 1 fields, both at loop order l; 
this appears to spoil the recurrence properties. The offending terms come from p = 0 (p = n) and 
p = 1 (p = n — 1). When p = 0 a tadpole multiplies ( d 2 + to 2 )( 1... n, x ), but it is only the tree 
level tadpole which multiplies ( d 2 + to 2 )( 1... n, x )i, and this term vanishes since by definition 
there is no tree level tadpole. The right hand side therefore contains correlation functions of n + 1 
fields but of order at most l — 1. 

When p = 1 there is a tree level two point function multiplying ( d 2 + m 2 )( 1... n — 1 ,x)i, 
but this becomes a delta function after integration by parts and we can move this contribution 
to the left hand side. The right hand side then contains correlation functions of n fields but of 
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order at most l — 1. We therefore have genuine recurrence relations in the number of fields and 
loop order. Following these details our result becomes 


I(n,l)( 1 ...n)i J Ax i(d 2 + m 2 ){ 1 ...n, x,x)i 


+ 2 S y m 


n — 2 l p 

EE/ 

v— n J 


Ax (l...p, x )i 


i(d 2 + m 2 


P= 2 ('=0 ‘ 
i-i 


p\(n — p)\ 

+ E E/ dx ( 1, x )i'i(dx + ra 2 )( 2 ... n, x )i_i 

1—>2—..n l' = 1 ^ 

Z-l - 

+ Ax (x)pi(d 2 + m 2 ){l.. .n,x)i-i>. 

Z / = 1 ^ 


(p+ 1.. .n, *);_(/ 


(17) 


At tree level there are no connected n-point functions for n < N and our result is empty in this 
sector. For n > N we have the simple relation 


(1.. .n)o = 


n — N 
2N - 4 


n — 2 ,, 

Sym ^ / da: (l...p, a: )o i(d 2 + m 2 )(p + 1... n, x)o- 

p=2 ■’ 


(18) 


In all of the expressions above we may insert under each integral a resolution of the identity, 

1 = J Ay i(d 2 + m 2 )G xy , (19) 

which can be used to extend the recurrence relations to scattering amplitudes. The tree level 
relations become, for example, 

^ 1 -" n ^ 0= 2N - 4 Sym E J dxd y [ i (dy+m 2 ){l...p,y)o\ p] [i(d 2 + m 2 ){p+l... n, x ) 0 ] . 

( 20 ) 

Multiplying by a product of operators i(d 2 +m 2 ) for each variable and taking the Fourier transform 
gives, on the left hand side, the expression for the n-particle scattering amplitude. For example, 
setting N = 3, n = 4 and l = 0 we recover the result for the four particle scattering amplitude in 
cf) 3 theory from the previous section. 


4 Quantum electrodynamics 


Our QED Lagrangian in the gauge — 0 is 

£ = iA M (9 2 77 MI/ ) A„ + ip(i-y.d — m)ip — etp^ipAp (21) 

which is obtained from the standard action by including the gauge fixing term 


<5 0 [WM„ + i6A] 

with BRST transformations 8qA ^ = <9 M c, SqA^ = 0, 5b = A, <5A = 0. The ghost integrations give 
an irrelevant factor to the partition function while the integral over A is Gaussian and removes 
the piece proportional to from the term quadratic in the gauge field. As before we will 

calculate the expectation value of the free action in the presence of sources, 


J T)(A, ip, tp) p J Ax ^A )1 (d 2 ri ll ' l ')A v +i/j(i'y.d - m)ip 


iS/h+iS src /h 
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where the source action is 


da; + %l) a r] a + 

for vector J and Grassmann valued (independent) spinors r/ and fj. Introducing the parameter £ 
the expectation may be written 




6 5 5 

+ih— 
oJ orj or] 


Z[J,r],ri\ 


d_ 

ac 


fv(A, 

C=o J 




+ iSint / b-^-iSsrc / b (*^2) 


Written in terms of functional derivatives the free action acts on W, the generator of connected 
graphs, as 


J da; (-iS 2 77 MI ') a 


SW 8W h 
8Jr(x) SJ"(x) + 2 
SW 


J da: (-id 2 ^):, 


5 2 W 


8Jr-^x)5J v (x) 


+ h [dx T—rri-1-dx - im)-™ + h [dx (-'y.dx - im) a 
J 5r(xy 5ri b (x) J 


S 2 W 


(23) 


5t] a {x)5r] b {x) 


The processes represented by these four terms are illustrated below, respectively. 




The first order derivatives join two diagrams by joining one external line from each to form a 
single internal line. Proving results beyond tree level has been one of the challenges of research 
in MHV and related areas. The terms with second order derivatives, above, deal with the loops 
in our approach. They act on a single diagram, joining two external lines together to form a 
loop, or cutting a loop open which decreases an amplitude’s loop number by one but adds two 
additional external legs, thus relating n and n + 2 particle scattering of one lower loop order. 

As we now have three different fields our equations will become more cluttered. The results 
will be clearer if we first calculate the recurrence relations between purely bosonic and fermionic 
amplitudes separately, by calculating the expectation value of only the bosonic or fermionic action 
with only bosonic of fermionic sources. The ( trick then counts only bosonic or fermionic lines. 
Since QED Feynman diagrams obey o and the additional constraint 

rif + 2 If = 2(nb + ‘21b) 

we can identify Jf, and If in terms of the number of external lines and the loop order, 

lb — g n f + ^ — 

If = n b + ^rif +2(1 - 1), ( 24 ) 

V = Tib + Tlf + 2(1 — 1). 
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This means that the effect of £ is still to multiply each < n >i by a function of only n and l. 
Expanding in the number of fields and the loop order, evaluating the first two terms in 12,11 . the 
recurrence relations for photon scattering are 

(n + Ib(n, 0)(A m (xi) ... A,j, n (xn))i = X J<Ax {-idlr] >J ' v )(A IJ . 1 (xi)... A fn (i„), A p (x), A„(x))i-i 

p=n 

1 ^ —id 2 71^ 

+ q Sym ^ ] (^mi O^i) • • • ( x p)i ,( — ^.(A Pp+1 (x p - |_i)... A Pri (x n ), A v ,(x))i_y. 

z p=Q P\n P). 

i '=o 

(25) 


The generating functional for purely fermionic amplitudes is 


Z[i1,V\ = exp ^S[ip,ip] + J dx rjip + ipri ). 


With this definition it may be checked that 
,2 n 8 5 8 


-Z[y,y] 


8r/ ai (xi) ''' 8r/ a » (x n ) <5p bl (pi) ''' 5r/ b ™ (y n ) 

- Jv(*p, (xi) ... ja n (x n )i, bl (wi)... A n (yn)ei sl ^ ]+ * I dx **+*> 

so that the logarithm of Z may be expanded as 

-t OO | /• oo 

= tXX TTjTJ / ■■■’n bl V an ■■■V ai ■■■ipbjlti 1 

n n =0 V V Z=0 


where the fermionic connected correlation functions have been abbreviated to 


(lp ai ■ ■■i>b n )l ■= (0 I T[lp ai (xi) . . ■'pa n (x n )lpb 1 {yi) ■ ..1pb n {yn)} | 0)i, 


and the time ordering observes the fermi statistics of the fields. The ( trick here multiplies the 
2n-field /-loop correlation function by — 2n — //(2n, l). The only changes to the bosonic case when 
evaluating the functional derivatives is to keep careful track of the minus signs. We find 


■h 


(2 n + If(2n, l)){ip ai ■ ..ipb n )i = - / Ax (-~/.d x - im) a (ip ai ... 'ipa n i’a.ip b ipb 1 ■ ■ ■ ipb n )i-i 

{-"T-dx - im) ab 


p=n 

l'=l 


L =L /• 

ASymXX y d:c (V^i ■ ■ • ‘>la p _ 1 Mb 1 ■■■ipi, p )i 


p =i 
i'=0 


K n 


-kPa p ■■■■4 , a n lpb1pb p+ 1 ...tpb n )l-l> 


(26) 


where K np = (—)d' 1_1 h p_1 )p!(p — l)!(n — p + l!)(n — p)! and ‘ASym’ is the instruction to anti 
symmetrise over the possible external field variables, with the ordering on the LHS given sign 
+1. 

Giving the full set of equations for all QED correlation functions is a simple extension, though 
we will have to reduce our notation further to present the result. Including all the sources in the 
partition function, we expand W as 

1 00 r 

W = j XX / R( n , rn)J^ 11 ... J^ n ff n ... rj bl Tf m ... r/ ai XX( 1 ■ ■ ■ n\l... m\l ... m)ih l 

n,m =0 1=0 
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where R(n,m ) = i n /(n!(m!) 2 ) and the correlation functions are 

(1 ■ • -n|l... m\l... m)i := (0 |T[^ M1 ... ■ ■■tpa m '^b 1 ■ ■ -ipb m \| 0 )z. (27) 

Equating by order in h and the number of fields gives 

(E + /)(1.. . n|l.. . m|l. ,.m)i = ^ J dx dl)( 1.. .n, {x,fi}, {x,v}\l.. .m\l .. 

- Jdx (—y.d x - im) ab ( 1... n\l.. .m, {a, *}|{6, x}, 1... m)i -1 


+Perm ^ dx(l. ..n, {x,n}\l . ...m' );/—-- 

/ n / ,/ n J ^Wn'm' 

n —0,m — l,i —0 

{n ... n, {x, v}\m! + 1... m|m , + 1... 

n,7n,l „ tri \ab 

—Perm ^ dx (1... n\l ... m! — 1, {a, ®}|1... m!)y —-—— 

/a / n J Vn'm' 

n' =0,m' = 1,1' =0 

(n + 1... n|m/ ... m\{b, x}m' + 1... 

where E + I = n + 2m + + If, Q n 'm' = (—) mm -m n'!(n — n')!(m'!) 2 (m — m')! 2 , V^/m/ = 

K mm /n'!(n — n')\ and ‘Perm’ is the instruction to symmetrise over the bosonic variables and anti 
symmetrise over the fermionic variables. 


5 Conclusions 

In this paper we have used the cutting and sewing of Feynman diagrams to describe the origins 
of recurrence relations in the correlation functions and scattering amplitudes they underlie. We 
have derived recurrence relations for the correlation functions in scalar held theory and quantum 
electrodynamics which hold to all loop orders. Inserting a resolution of unity into the relations, 
as in the scalar held case m . generates recurrence relations between scattering amplitudes. 

The Schwinger-Dyson equations can be used to give an alternative derivation of our results 
as we will now describe. The QED Schwinger-Dyson equations for correlation functions of the 
photon held are 


i + SS ^,Tp] + jA^/KHU^/K _ Q 

Taking a derivative with respect to the source gives 

Tr ( v )5{0)Z[J]+ l - (^S 0 [A] + Si^A,^,^] - ihj^-^e iS/h+i * A » /h = G. (28) 

The hrst term is divergent and should be regulated, but we will see in a moment that it cancels 
another divergence coming from the free piece of Z[J\. Using the £ trick from earlier we know 
that the third term in counts the number of vertices in a correlation function. By expanding 
W in powers of J it is easily checked that the fourth term in (HU counts the number of external 
lines. We will evaluate the second term in terms of functional derivatives as we did before, but 
this time we will isolate the free part of the generating functional, defining W, 

W[J ] = ~ J + w[j}. 
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(29) 






Evaluating the free action in terms of functional derivatives we find that the divergence cancels 
against a similar term from the free field partition function and the remaining terms depend only 
on W. The result is 


■/ 


hi dx (—*9^'“') 


8 2 W 


= J dx J M ( x) 


5 J v (x)8 J v (x) 8J^{x 

8W 1 


2 SW 

(—tOxV )- 


5J^(x) h 


8J v {x ) 

OO p - n oo 

E / h jn Y, h l NV{n,l){l ...n)i 

^ —n J ' i — n 


(30) 


which is equivalent to our earlier result for the scattering of photons, equation (1551 . The same 
reasoning with the scalar field directly leads to a tidier version of im in terms of W. 

Our method also gives us an elegant interpretation of the Schwinger-Dyson equations in 
perturbation theory. Consider the scalar held Schwinger-Dyson equation, 


i 

h 



&So[0] + + Jf x \\ f J<t>/h _ q 

84>(x) 8(j>{x) ) 


Taking a derivative with respect to J(x) gives 


S(0)Z[J] + ^j D<^2So[</>] + NSi n t[<f>] - ih J J (x) e iS/h+i 1 = 0. (31) 

As before the divergent first term cancels a divergence in the free part of Z and we will neglect 
it. The first term in large brackets counts minus twice the number of internal and external lines, 
—2(n + I), in a diagram. The second term counts N times the number of vertices, NV. The 
third term counts the number of external lines. Therefore the non-perturbative Schwinger-Dyson 
equation is equivalent to 

—2(n + I) + NV + n = 0, 

in perturbation theory, which is just the Euler equation. 

The relations we found are of the same, very basic, structure found in the MHV rules of Yang- 
Mills theory - diagrams for n particle scattering may be separated into diagrams for lower order 
scattering by severing an internal (non loop) line, or built from n + 2 particle diagrams of one 
lower loop order by attaching two external legs together to form a loop. Although our approach 
does not give a reason why amplitudes with particular helicities should have special properties, 
it does suggest that the recurrence relations hold in a much wider context. A issue with our 
method, however, is that in theories with interactions of multiple orders, such as Yang-Mills and 
light cone string Held theory, the number of vertices and internal lines is no longer unique in 
the diagrams of a correlation function so the effect of the parameter £ no longer has the simple 
interpretation of im In principal the extension is very simple; the expectation value of the whole 
action gives a counting which is unique in the above sense but requires evaluating the interaction 
terms as functional derivatives which would give much more complex results. After posting this 
paper we became aware of uni, where similar arguments with the Schwinger Dyson equations 
were used to investigate precisely this issue in Yang-Mills theory. The recursion relations between 
scattering amplitudes given in that paper are similar to our own but contain terms cubic as well 
as quadratic in the amplitudes which originate in the inclusion of interaction terms. However, 
simple diagrammatic examples such as that used in chapter 2 can be constructed for theories with 
interactions of multiple orders. Due to this and the MHV results in the literature we suspect 
that the effects of the interactions on the recursion relations may be largely degenerate. We will 
address this in a future publication; in particular we will consider the scattering of strings as 
described by string field theory. 
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